Charge Induced Vortex Lattice Instability 
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It has been predicted that superconducting vortices should be electrically charged 
and that this effect is particularly enhanced for high temperature superconductorsP^ 
Hall effect 3 and nuclear magnetic resonance (NMR) experiments 4 suggest the existence 
of charge accumulation in the vortex core, but the effects are small and the interpre- 
tation controversial. Here we show that the Abrikosov vortex lattice, characteristic 
of the mixed state of superconductors, will become unstable at sufficiently high mag- 
netic field if there is charge trapped on the vortex core. Our NMR measurements of 
the magnetic fields generated by vortices in Bi2Sr2CaCu2 0s+ y single crystals 5 provide 
evidence for an electrostatically driven vortex lattice reconstruction with the magni- 
tude of charge on each vortex pancake of ~ 2xl0 _3 e, depending on doping, in line with 
theoretical estimatesP^ 



The behavior of quantized vortices in high tempera- 
ture superconductors (HTS) affect many of their appli- 
cations ranging from high current carrying wires^ to elec- 
tronic devices P Effective performance requires that the 
vortices be stable since their dynamics lead to dissipation 
originating from the electronic excitations within the vor- 
tex cores. Effort to understand the interactions between 
vortices has been focused on the circulating supercur- 
rents that are characteristic of the vortex and generate 
local magnetic fields. 9 However, if the vortices are electri- 
cally charged there will also be an electrostatic interac- 
tion, possibly relevant to vortex stability. We have per- 
formed model calculations showing that if vortices trap 
electrical charge, then the Abrikosov vortex lattice be- 
comes unstable in sufficiently high magnetic fields where 
the magnetic interactions diminish and the electrostat- 
ics dominate. Our NMR experiments confirm the exis- 
tence of a magnetic field induced instability in the vor- 
tex structure of the highly anisotropic superconductor 
Bi 2 Sr 2 CaCu 2 8+2y (Bi2212). 

A strong magnetic field normal to the conduction 
plane, H\\c, of Bi2212 penetrates in the form of quan- 
tized flux bundles, called pancake vortices. 5 They inter- 
act magnetically, with a repulsive interaction for pan- 
cakes within a plane, but attractive if they lie in different 
planes. The attractive force is responsible for aligning the 
vortex cores, one above the other, thereby forming flux 
lines as indicated in Fig. 1. Perfect alignment is called 
the Abrikosov configuration. However, the average at- 
tractive force decreases^ with increasing field as the vor- 
tex density increases and as the magnetic field modula- 
tion from nearby vortex supercurrents cancel each other 
more effectively. In contrast, for charged vortex cores, the 




FIG. 1: Vortex pancakes with competing interactions. 

Attractive- magnetic and repulsive-electrostatic forces com- 
pete and change dominance as the magnetic field is increased. 
The magnetic interactions between vortices in different planes 
is attractive preferring a line of vortex pancakes while the 
Coulomb interaction is repulsive. 



electrostatic force between pancakes on different planes 
increases at short range. A complete picture of the bal- 
ance of these forces requires that the nearby stacks of 
vortices be taken into account. Clem calculated the cost 
in magnetic energy for displacing an entire plane of vor- 
tices, relative to those in all other planes P We follow 
this approach and compare the magnetic energy with the 
change in electrostatic energy when there is a charge on 
each vortex core. Such displacements can become more 
favorable for charged vortices, decreasing the total en- 
ergy as is illustrated in our calculations shown in Fig. 2. 
At sufficiently high magnetic field the electrostatic and 
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FIG. 2: Energy phase diagram. Two separate regions in 
the charge-field plane correspond to dominance of charge or 
magnetic forces. Below the red curve the magnetic energy is 
larger and above the red curve the electrostatic energy domi- 
nates. The dashed lines are chosen to intersect this phase di- 
agram at the magnetic field where the NMR second moments 
have a minimum, Fig. 3. [Inset] At each magnetic field, for a 
charge of 2.1xl0 _3 e, the change in total energy is shown for 
a rigid displacement of a layer of vortices in the ab plane to 
various positions ranging over several vortex unit cells of the 
Abrikosov configuration. Softening of this energy landscape 
is evident at 10 T. 



magnetic interactions balance and there is an instability. 
The effect is evident in the inset of this figure at a field 
of 10 T for a charge of 2.1xl0 _3 e, where the energy, as 
a function of vortex layer displacement in the ab plane, 
becomes fiat. 

The existence of charge trapped on the vort ex cor e of a 
superconductor has been discussed extensivel y 1 ^^K lU as 
a consequence of the Bernoulli potential for the electro- 
static fields, and as early as 1961 by Fritz London. 16 Even 
for conventional superconductors the electrodynamics of 
vortices lead to a depletion of charge in the vortex core in 
order to balance inertial and Lorentz forces.^ Addition- 
ally, if there is particle-hole asymmetry, it is favorable 
for charges in the core states to move to lower energy 
in the outer superconducting region in order to main- 
tain constant chemical potential!^ These di fferen t mech- 
anisms produce a vortex charge per pancake^^ of order 
Q ~ e(A/ei?) 2 , a reduction from the electronic charge, e, 
by the ratio of the superconducting gap to the Fermi en- 
ergy squared. This can be of order 10 _6 e for low temper- 
ature superconductors but more favorably, ~ 10 _3 e for 
HTS, which have a larger gap. Vortex charging ha s been 
studied theoretically for d-wave superconductor^ 6 -* 17 * 18 * 
coupled to antiferromagnetic (AF) order in the vortex 
core, and chiral triplet superconductors.^ The experi- 
mental situation is not so well-established. A sign change 
in the Hall effect 3 has been interpreted 1 as evidence for 
vortex charging. However, sign changes have been re- 
cently ascribed to Fermi surface effects Comparison 
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FIG. 3: Softness of the charged vortex structure. The 

dashed curves are the root-mean-square energy cost for vortex 
displacements from the Abrikosov configuration perpendicu- 
lar to the applied field for the effective vortex charges of 2.1 
and 2.5 x 10~ 3 e. These curves (left axis) are shown super- 
posed on our NMR data, to demonstrate that the softening 
of the vortex- vortex interactions for these charges, occurs at 
the magnetic fields where the NMR second moment is a min- 
imum. The second moments of the 17 NMR spectra (right 
axis) are shown for three crystals independently processed by 
17 exchange to give T c — 75 K (blue squares) and T c = 82 
K (red circles, open and closed), at T = 4.2 K (error bars 
represent the background noise). The data are normalized to 
the theoretical value 10 for a Abrikosov vortex configuration 
at low magnetic field, ao = 0.0037100 /A 4 . An independently 
determined 5 Knight shift distribution that contributes 1.25 
kHz/T to the linewidth, has been subtracted. At H — 8 T 
this corresponds to 0.44 of the normalized second moment. 



of the nuclear quadrupole resonance frequency and the 
satellite splittings in NMR has also been interpreted as 
evidence for vortex charging.^ But the magnitude and 
sign of the charge, and the temperature dependence are 
in disagreement with BCS theor^U^ a s well as d-wave 
theoryP 

We model the competing magnetic and electrostatic 
interactions by calculating the total energy per vortex 
for a rigid displacement of a single layer of pancake 
vortices in the ab plane relative to the energy of the 
Abrikosov configuration. Josephson coupling between 
planes is suppressed at high magnetic field by temper- 
ature fluctuations^ an d can be ignored. From Clem's 
calculation,^ the magnetic coupling energy per vortex is 
given by, 



E m (r) = (H/M Uo(s)[e is - r ~ 1] (1) 



where, 



U {q) = -^sl^\% 2 + \- 2 ). (2) 



Here H is the applied field, normal to the conducting 
planes; g is a reciprocal lattice vector of the Abrikosov 
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FIG. 4: Electrostatically driven phase change of the 
vortex lattice. The red symbols, both filled and open, rep- 
resent vortex pancakes in the Abrikosov configuration. At 
magnetic fields above the instability three inequivalent lat- 
tices appear (red, blue, and green). The aspect ratio in 
the bottom portion of the figure matches the geometry of 
Bi2Sr2CaCu208+y at H = 10 T. As the magnetic field in- 
creases the Coulomb potential becomes greater than the mag- 
netic inter-layer coupling, making it favorable for pancake 
vortices to displace as far as possible from near neighbors 
in adjacent planes. 



vortex lattice in g-space; 0o the flux quantum; A = 220 
nm is the in-plane penetration depth 5 ; the bilayer spacing 
is s = 15 A; and r the displacement of the plane of pan- 
cakes restricted to the vortex unit cell of the Abrikosov 
configuration (r = 0) . The sums converge when taken to 
seven times the vortex spacing. For the electrostatic en- 
ergy we calculate the Coulomb potential of a periodic lat- 
tice and determine the energy required for the displace- 
ment defined above. The Coulomb potential for point 
charges is summed in direct space assuming a cut-off of 
35 times the inter-vortex spacing. To maintain charge 
neutrality, there is a positive background charge spread 
throughout the unit cell, but when considering the energy 
difference of a displacement, its effect is canceled. The 
calculation of the energy landscape for the displacement 
of a single plane is displayed for various magnetic fields 
in the inset to Fig. 2. At low magnetic fields the lowest 
energy corresponds to the Abrikosov configuration where 
all vortex cores line up, one above the other. However, 
in high magnetic fields, the most favorable location is for 
vortices in the displaced plane to line up at the mini- 
mum magnetic field position equidistant from vortices in 
the plane above it (or below it), as sketched in Fig. 4. 
It is remarkable that at an intermediate field the energy 
landscape is essentially flat. This condition can be ex- 
pressed quantitatively in terms of the root-mean-squared 
energy averaged over the unit cell, shown in Fig. 3. The 
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FIG. 5: Magnetic field dependence of the second mo- 
ment for the vortex local-field distribution. The black 
curve is calculated from Ginzburg Landau theory for the 
Abrikosov configuration.^ The green curve corresponds to 
a tripled vortex unit cell which in our model is an electrostat- 
ically driven reconstruction of the vortex lattice obtained by 
changing 0o to 0o /3 as suggested by BrandtP^ and depicted in 
Fig. 4. There is a factor of 2 to 3 between the second moments 
of the two configurations at any given field. Such a vortex lat- 
tice reconstruction leads to a decrease in the second moment 
similar to our data, Fig. 3. [Inset] The calculated NMR spectra 
at 10 T for both configurations of vortex pancakes are plotted 
with respect to their first moments. 



minimum in this rms-energy as a function of magnetic 
field defines the transition between magnetic and electro- 
static order for a given vortex charge. The calculations 
are performed for different charges giving the phase di- 
agram in Fig. 2. This calculated charge is really an ef- 
fective charge which includes the Thomas-Fermi screen- 
ing within a plane responsible for the suppression factor, 
(A/e^) 2 , as noted above. But it does not allow for di- 
electric screening that will increase the charge by 
where the dielectric constant, ft, is a materials parameter 
~ 5) 14 l 21 l Since Bi2212 is a bilayer system, the charge per 
conduction plane in a vortex is larger than the effective 
charge by ~ 1. 

The experimental NMR spectrum directly reveals the 
distribution of magnetic fields generated by vortices and 
their spatial correlations. In the superconducting state 
the spectrum is simply proportional to th e pro bability 
distribution of the local fields. It is commoiP^ to com- 
pare experimental spectra with that expected for a regu- 
lar array of vortex lines. In this case the field distribution 
has a characteristic asymmetry and singularities associ- 
ated with the minimum local field, the saddle point in 
the field distribution, and a maximum field at the vortex 
core, inset to Fig. 5. The second moment of the spec- 
trum, closely related to the square of the linewidth, is 
a direct measure of the inhomogeneity of the fields from 
vortex supercurrents and can be interpreted in terms of 
the superconducting penetration length, A, and the ge- 
ometry of the vortex array. If the vortex pancakes are 
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disordered within the conduction plane, the distribution 
is broadened. In contrast, disorder or reconfiguration 
of vortices from plane to plane decreases the linewidthP^ 
From our NMR experiments we find an unusually abrupt 
decrease in the linewidth with increasing applied mag- 
netic field for 17 NMR measured at T = 4.2 K at the 
0(1) oxygen site in the copper-oxygen conducting plane 
of Bi2Sr2CaCu208+ 2/ . The magnetic field dependence of 
the second moment of the spectrum is shown in Fig. 3. 
The magnitude of the second moment at lower fields, 
H < 5T, is of order that expected from an Abrikosov con- 
figuration. We interpret the linewidth collapse with in- 
creasing field as a transformation to a new, field-induced, 
vortex state. 

The abrupt decrease of the NMR linewidth with in- 
creasing magnetic field for several over-doped crystals 
(T c = 75 K, 82 K), Fig. 3, provides consistent experimen- 
tal evidence that charged vortices can be a mechanism 
for vortex structure instability. Additionally, our results 
indicate that this instability moves to higher magnetic 
fields with increased doping, corresponding to a smaller 
vortex charge as suggested by calculations for a d-wave 
superconductor. 6 In related work, Gurevich found that 
if a single vortex line is charged, it becomes unstable to 
chiral distortion.^ 

In our model the Abrikosov configuration is trans- 
formed to an A-B-C-type of tripling of the vortex unit 
cell along the c-axis as sketched in Fig. 4. At high 
magnetic fields a pancake vortex seeks out a position 
that is anti-correlated with vortices in planes above and 
below it. From the minimum NMR linewidth in Fig. 3 
we find the magnetic field for this transition. Then 
according to the phase diagram in Fig. 2, the effective 
charges per pancake are 2.1 and 2.5 x 10 _3 e for the 
crystals T c = 75 K and T c = 82 K, respectively. This 
new ordering maintains the intra-plane vortex-vortex 
correlations but the spacing of vortex cores on adjacent 
planes is decreased by ~ 1 / y/3. This leads to a decrease 
in the distribution of local magnetic fields as pointed 
out by BrandtpSI and a corresponding decrease in the 
second moment of the NMR spectrum as we show in 
Fig. 5. Raw spectra and further experimental details are 
presented in the Supplementary Materials section. On 
close examination the NMR spectra at low field have two 
distinct components which merge with increasing field. 
The composite spectra have an overall shape rather 
similar to the superposition of the calculated spectra 
shown in the inset to Fig. 5. We infer that it is likely 
that the vortex rearrangement proceeds continuously 
over a range of magnetic field and completes at the 
field for which the spectrum width is a minimum. The 



minimum is a balance between the abruptly decreasing 
vortex contribution to the local field distribution and 
an increasing paramagnetic contribution with increasing 
applied field that we have followed up to as high as 30 
T, Fig. 3. We identify this high field behavior with the 
magnetization induced by bound states in the vortex 

core J23] 



I. METHODS 

Our samples were prepared by exchanging in 1 bar of 
flowing 17 (70-90% enriched) at 600° C for 48 h followed 
by annealing for 150 h at 450° C. Our 82 K data in 
Fig. 3 are from two separate samples from the University 
of Tokyo. The closed circles correspond to three single 
crystals of Bi-2212, stacked along the c-axis, with a total 
mass of 20 mg. The sample represented by the open 
circles has a similar total mass, though the crystals were 
cleaved to improve signal-to- noise ratio. The squares are 
data from a previous experiment^ with crystals, T c = 
75 K, provided by Argonne National Laboratory with 
similar processing conditions. 

NMR measurements were performed in magnetic fields 
from 4-30 T at a temperature of 4.2 K at Northwestern 
University and the National High Magnetic Field Lab- 
oratory in Tallahassee, Florida. We have shown that at 
this temperature the vortex structure is frozen at all mag- 
netic fields P Spectra were taken using a frequency sweep 
technique with a typical ir/2 pulse of 4 /is. Before tak- 
ing data at any particular magnetic field, the sample was 
heated above T c for 10 min and then quench-cooled to 
4.2 K. This method produces narrow spectra, with mini- 
mal flux pinning, and was found to give consistent results 
through multiple field cycles. In general flux pinning is 
weaker at higher magnetic fields, and if it exists can be 
identified using field-sweep NMR. 



II. ACKNOWLEDGEMENTS 

We thank E.H. Brandt, L.N. Bulaevskii, C.A. Collett, 
M. Eschrig, M.R. Eskildsen, W.J. Gannon, A. Gurevich, 
A.E. Koshelev, Jia Li, V.F. Mitrovic, J. Pollanen, and 
J. A. Sauls for helpful discussions. This work was sup- 
ported by the Department of Energy, contract DE-FG02- 
05ER46248 and the National High Magnetic Field Lab- 
oratory, the National Science Foundation, and the State 
of Florida. 



1 Khomskii, D. I. & Freimuth, A. Charged vortices in high 2 Blatter, G., Feigel'man, M., Geshkenbein, V., Larkin, A., 
temperature superconductors. Phys. Rev. Lett. 75, 1384- & van Otterlo, A. Electrostatics of vortices in type-II su- 

1386 (1995). per conductors. Phys. Rev. Lett. 77, 566-569 (1996). 



5 



3 Hagen,S. J., Lobb, C. J., Greene, R. L., & Eddy, M. Flux- 
flow Hall effect in superconducing TbE^CaC^Os films. 
Phys. Rev. B. 43, 6246-6248 (1991). 

4 Kumagai, K.-L, Nozaki, K. & Matsuda, Y. Charged 
vortices in high-temperature superconductors probed by 
NMR. Phys. Rev. B 63, 144502 (2001). 

5 Chen, B. et al. Two-dimensional vortices in superconduc- 
tors. Nature Physics 3, 239-242 (2007). 

6 Knapp, D., Kallin, C, Ghosal, A. & Mansour, S. Antifer- 
romagnetism and charged vortices in high-T c superconduc- 
tors. Phys. Rev. B 71, 064504 (2005). 

7 Shen, T. et al. Development of high critical current den- 
sity in multifilamentary roundwire Bi2 Sr2CaCu2 0s+(5 by 
strong overdoping. Appl. Phys. Lett. 95, 152516 (2009). 

8 Martinis, J. M., Ansmann, M. & Aumentado, J. Energy 
decay in superconducting Josephson-j unction qubits from 
nonequilibrium quasiparticle excitations. Phys. Rev. Lett. 
103, 097002 (2009). 

9 Clem, J. R. Two-dimensional vortices in a stack of thin 
superconducting films: A model for high-temperature su- 
perconducting multilayers. Phys. Rev. B 43, 7837-7846 
(1991). 

10 Brandt, E. H. Magnetic-field variance in layered supercon- 
ductors. Phys. Rev. Lett. 66, 3213-3216 (1991). 

11 LeBlanc, M. Vortices in high temperature superconduc- 
tors: role of the centrifugal and Lorentz forces. Supercond. 
Sci. Technol. 10, 919-931 (1997). 

12 Kolacek, J., Lipavsky, P. & Brandt, E. H. Charge profile 
in vortices. Phys. Rev. Lett. 86, 312-315 (2001). 

13 Simanek, E. Charged vortices and copper nuclear 
quadrupole resonance in the cuprates. Phys. Rev. B 65, 
184524 (2002). 

14 Lipavsky, P., Kolacek, J., Morawetz, K. & Brandt, E. H. 
Electrostatic potential in a superconductor. Phys. Rev. B 
65, 144511 (2002). 

15 Eschrig, M. & Sauls, J. A. Charge dynamics of vortex 
cores in layered chiral triplet superconductors. New J. Phys 
11, 0750 09 (2009). URL |http://iopscience.iop.org/1367-| 
2630/11/7/075009/, 

1D London, F. Superfluids (Dover, New York, 1961). 

17 Chen, Y., Wang, Z. D., Zhu, J.-X. & Ting, C. S. Vortex 
charges in high-temperature superconductors. Phys. Rev. 
Lett. 89, 217001 (2002). 

18 Zhao, H.-W., Zha, G.-Q., Zhou, S.-P. & Peeters, F. M. 
Long-range Coulomb repulsion effect on a charged vor- 
tex in high-temperature superconductors with competing 
d-wave and antiferromagnetic orders. Phys. Rev. B 78, 
064505 (2008). 

19 LeBoeuf, D. et al. Electron pockets in the Fermi surface of 
hole-doped high-T c superconductors. Nature 450, 533-536 
(2007). 

20 Glazman, L. I. & Koshelev, A. E. Thermal fluctuations 
and phase transitions in the vortex state of a layered su- 
perconductor. Phys. Rev. B 43, 2835-2843 (1991). 

21 Maki, M. et al. Dielectric properties and carrier dynam- 
ics in Bi2Sr2YCu208+<5- Journal of the Physical Society of 
Japan 76, 044711 (2007). 

22 Brandt, E. H. Precision Ginzburg-Landau solution of ideal 
vortex lattices for any induction and symmetry. Phys. Rev. 
Lett. 78, 2208-2211 (1997). 

23 Mitrovic, V. F. et al. Spatially resolved electronic structure 
inside and outside the vortex cores of a high-temperature 
superconductor. Nature 413, 501-504 (2001). 

24 Gurevich, A. Helical instability of charged vortices in lay- 



ered superconductors. Phys. Rev. B 81, 100505 (2010). 



III. AUTHOR CONTRIBUTIONS 

Experiments were carried out by A.M.M., S.O., S.M., 
W.P.H., A.P.R. and P.L.K. Samples were provided by 
K.F., M.L, and S.U. 



IV. AUTHOR INFORMATION 

Correspondence and requests for materials should be 
addressed to w-halperin@northwestern.edu. 

V. COMPETING FINANCIAL INTERESTS 

The authors declare no competing financial interests 

VI. SUPPLEMENTARY MATERIALS 
A. Bi-2212 17 NMR 




14.95 



14.90 14.85 
Magnetic Field (T) 



14.80 



FIG. 6: The 17 NMR spectrum of the T c = 82 K 
sample O(l) site at a frequency of 86 MHz and tem- 
perature 4.2 K. The central transition (±1/2 <H> =Fl/2) 
is shown along with the first pair of quadrupolar satellites, 
(±3/2 <H> ±1/2), with a quadrupolar splitting of Aisq = 63 
kHz. For the orientation of our sample, H || c-axis, the central 
transition has no quadrupolar broadening. 

With multiple stoichiometric locations and quadrupo- 
lar resonances, the full oxygen NMR sectrum of Bi-2212 
is complex. The main peak in Fig. 6 is the central tran- 
sition, (±1/2 o =Fl/2), of the 17 spectrum of Bi-2212 
from the 0(1) site. This site is understood to be the oxy- 
gen resonance from the superconducting Cu02 planes. 
The 0(2) resonance from the oxygen in the insulating 
SrO planes is completely saturated in our measurements 
at 4.2 K and is not discernible in Fig. 6. The oxygen 
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in the BiO plane is unobservable due to structural disor- 
der in this plane. The 0(1) spectrum, in addition to 
the central transition, is split into two different pairs 
of quadrupolar satellites for the (±3/2 <o> ±1/2) and 
(±5/2 o ±3/2) transitions, the latter are just outside 
of the frequency range presented in Fig. 6, see Ref 5 for 
the full spectrum. 
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FIG. 7: Field evolution of the central O(l) NMR spec- 
tra for crystals with T c = 75 and 82 K. Each spectrum 
was recorded after the sample was field cooled through the 
transition temperature to avoid flux pinning as described in 
the Methods Section. The lineshape of the spectra at these 
low fields appears to have two components that merge near 
H = 4.5 and 6 T, respectively, for these two crystals. 

As the sample is cooled from the normal state, through 
the superconducting transition, and finally to the vortex 
solid phase, the NMR spectrum shows a clear temper- 
ature dependence described in Ref. 5. At relatively 
high magnetic fields and at temperatures well below 
T c , vortices form a solid phase; the spectrum broadens 
and becomes asymmetric as would be expected for the 



Abrikosov vortex lattice, Ref. 5. 



B. Multicomponent Spectra 

The broad spectra at low fields, for which the second 
moments are shown in Fig. 3, appear to be composite 
spectra with two principal components. These spectra 
are displayed in Fig. 7 for the two crystals with, (a) T c 
= 75 and, (b) 82 K. The spectra narrow as the lower 
frequency contribution decreases in intensity with pro- 
gressively increasing magnetic field. The sample with 
T c = 75 K has two components which are somewhat nar- 
rower and spread further apart, as compared with the 
crystal with T c = 82 K. We note that the calculated field 
distribution for two vortex states: the Abrikosov config- 
uration of aligned vortices and that for the anti-aligned 
vortex pancakes, have a shift in frequency in addition 
to a different magnetic field distribution, as displayed 
in the inset to Fig. 5. Their superposition is similar to 
the measured spectra in Fig. 7, suggesting that the in- 
stability observed is a crossover from an inhomogeneous 
distribution of such vortex states to a more stable distri- 
bution at high magnetic fields. This complex behavior 
is not captured by our model for competition between 
electrostatic and magnetic interactions between vortex 
pancakes in different conduction planes. Nonetheless, it 
is qualitatively consistent with the general picture which 
we propose of a field-induced instability attributable to 
vortex charge. 



